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A model to describe the mechanism of conformational dynamics in secondary protein 
based on matter interactions is proposed. The approach deploys the lagrangian method 
by imposing certain symmetry breaking. The protein backbone is initially assumed to 
be nonlinear and represented by the Sine-Gordon equation, while the nonlinear external 
bosonic sources is represented by interaction. It is argued that the nonlinear source 
induces the folding pathway in a different way than the previous work with initially 
linear backbone. Also, the nonlinearity of protein backbone decreases the folding speed. 
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1. Introduction 

It is known that the nonlinear excitations play an important role conformational 
dynamics. For instance, the effective bending rigidity of a biopolymer chain could 
lead to a buckling inst ability Sorn e models have then been proposed to explain 
such protein transition I 2 I 3 | 4 | 5 | 6 | 7 | 8 | 9 ] 

In our previous work, the conformational dynamics of secondary protein can 
be modeled using 0'* interactions 10 11_ ^ j^g^g been shown that the model has 
reproduced the toy ad-hoc model based on the set of nonlinear Schrodinger (NLS) 
and nonlinear Klein-Gordon equation of motions (EDMs) in a more natural way 
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from first principle . In the model, the unfolding state of protein is initially assumed 
to be linear. On the other hand, the folding pathway is induced by the nonlinear 
sources like laser. Both protein conformation changes and the injected non-linear 
sources are represented by the bosonic lagrangian with an additional (j)^ interaction 
for the sources. It has been argued that the 'tension force' which enables the folded 
pathway can be reproduced naturally in the EOM. 

In this paper, we consider the nonlinear unfolding protein at the initial state, 
while the external sources remain nonlinear as previously done. This is important 
to investigate whether the folding mechanism and speed are influenced by the initial 
conformational state. Otherwise one cannot determine if the folded pathways are 
really induced and dominated by the nonlinear sources or not. 

The paper is organized as follows. First the model under consideration is pre- 
sented briefly, followed by the numerical analysis of the EOM induced by the model. 
Finally it is summarized by short discussion of the results and subsequent conclu- 
sion. 



2. The model 



This wor k is an extension of the previous model on protein folding using lagrangian 
approach contrast with the previous model which assumes the initial con- 

formational state is linear, now the protein is initially assumed to be nonlinear likes 
Sine-Gordon soliton. 
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(1) 



However, the sources injected into the backbone remain nonlinear and massless. 
Then, same as before the nonlinear sources are modeled by ijj'^ self-interaction. 

^. = ^M^(5^^) + ^(V'W- (2) 
The interaction term between both is described by, 

A„t = -A(0t0)(^t^) . (3) 

All of them provide the underlying model in the paper with total potential, 



^tot(^,0) 



1 




A 
4! 



^{i>Hf-K{cj,'^mH). (4) 



Now, throughout the paper let us assume that A^ is small enough, that is ap- 
proximately at the same order with A^ . In this case, the first term can be expanded 
in term of -v/A^, 



At. 



A 
4! 



^(V'W-A(<^t^)(^t^). (5) 
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up to the sec ond order accuracy. If = 0, the result coincides to the case in the ear- 
her wor Imposing namely local U(l) symmetry breaking to the total lagrangian 
makes the vacuum expectation value (VEV) of the fields yields the non-trivial so- 
lutions. In the preceding model, the 'tension force' which plays an important role 
to enable folded pathway has appeared naturally by concerning the minima of total 
potential in term of source field 



W = ('/')• (6) 



Beside of that, in term of conformation changes field the VEV is, 



6ml - 12A(V')2 

<*> = V V ■ p' 

It shows that the existence of Sine-Gordon potential makes the early stable ground 
state of conformational field turns out to be metastable. In other words, the non 
trivial VEV in Eq. (O constitutes new more stable ground state of the conforma- 
tional field. Transition between metastable into stable state breaks the symmetry of 
the vacuum spontaneously, while the conformational field should be nonlinear even 
though the external nonlinear source has not been instilled. Therefore the protein 
backbone should be in nonlinear form at the initial stage. 

The symmetry breaking at the same time shifts the mass term of <j) as follow, 

24A^ 

ml^ml = ml-^{4,f . (8) 

Nevertheless, the nonlinear source field is set being massless, since it represents a 
bunch of light source like laser. Thus, the broken symmetry of conformational field 
should not be considered to introduce its mass. 



3. The EOMs 

Having the total lagrangian at hand, one can derive the EOM using the Euler- 
Lagrange equations, 

where £tot = Cc + + ^int in Eqs. ([1]), ^ and ^ respectively. 

Substituting Eqs. ([1]), © and ([3]) into Eq. one immediately obtains a set of 
EOMs, 



''^1*1 '^^l*l.2A|«f-^si„^wUo, (10) 



2A|^||</.|^-:^|V'r = 0. (11) 
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Here the natural unit is restored to make the Ught velocity c and h reappear in 
the equations. The last terms in Eqs. (fTO|) and pT|) determine the non-linearity 
of backbone and source respectively. Also, the protein mass term is melted in the 
Sine-Gordon potential. One should put an attention in the second last term of Eq. 
pO)). i.e. k(j) with k ^ 2A{tp)'^. This actually induces the tension force which is 
responsible for the dynamics of conformational field and enabling the folded pathway 
as expected. 

Hence, solving both EOMs in Eqs. (fTO)) and ()11|) simultaneously would provide 
the contour of conformational changes in term of time and one-dimensional space 
components. The EOMs will be solv ed num erically using forward finite difference 
method as done in the previous work HillEl 



4. Numerical solution of EOMs 

Same as before, it is more convenient to replace -0 and (f> with u and w respectively 
and rewritten it in explicit discrete forms as follows. 



Mj+i = ^"i J - Uij^i -ret y h /.I'LW^jUi.j 

6 



<j] ' (12) 



■ I ~Z — : (l-j) 



for i = 2, 3, • • • , — 1 and j = 2, 3, • • • ,M — 1. Forward iterative procedure of the 
discrete EOMs can be performed if the two lowest time values are known. First, the 
value at ti is fixed by the following boundary conditions, 

■0(0, t) = ip{L, t)^0 and 0(0, t) = 0(L, t) = for < < < & , 

il^iXjO) — f{x) and (/)(x, 0) = p(x) for < x < L , ^^^-^ 

^^=,(.) and forO<.<L, 

with /(x), p{x), g(x) and q{x) are newly introduced auxiliary functions. Secondly, 
the values at can be determined using second order Taylor expansion, 

C^e^ f fi + l ~ '^fi + fi-l , r,A„2_f '^4' r3 



c^e^ ( p,+i - 2p, + p,^i 2 
■Wi,2 =P^- eg, + I ^ h 2A/, 



,3 , (16) 
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Fig. 1. The soliton propagations and conformational changes on the protein backbone inducing 
protein folding. The vertical axis in soliton evolution denotes time in second, while the horizontal 
axis denotes its amplitude. The conformational changes are on the (x, y, z) plane. The constants 
of the simulation are chosen as m = 0.08 eV = 1.42 X lO"''^ kg,L = 12 eV~^ = 2,364 nm,A = 
2.83 X 10-3, = 5 X 10-3, = 6 x 10-3,and h = c= 1. 

for i = 2,3, ■ ■ ■ , N — 1. S = Ax and e — At constitutes the side length between the 
discretized value. 

At the initial stage, suppose the nonlinear source and conformation fields have 
particular form of f{x) — 2sech(2a;) e*^^ and g{x) — arctan[exp(4a; — 10)], while 
g{x) = q{x) — for the sake of simplicity. Furthermore, the numerical solutions can 
be obtained by iterative procedure against Eqs. (|12p and using the results in 
Eqs. ((T5)) and (flB)) with the boundary conditions in Eq. (fT4|) . 

The procedure has been performed numerically and the results are given in Fig. 
[TJ The left figures in each box describe the propagation of nonlinear sources in 
protein backbone, while the right ones show how the protein is folded according 
to the time evolution. From the figure, it is clear that the protein backbone is 
infinitesimally bending at the initial stage before the nonlinear source injection. The 
bending constitutes the contribution of Sine-Gordon potential into the conformation 
field. However, this bending is too small to generate folding pathway, then the 
backbone still remains unfolded. 

The conformation changes which generate the folding pathway start appearing 
as the soliton starts propagating over the backbone. The result is surprisingly, even 
slightly, different with the earlier work. The folding processes are slower than the 
linear conformation case It might be considered as an effect of the nonlinear 
conformational field. One may conclude here that the effect is destructive against 
the nonlinearity of nonlinear sources. It can also be recognized from Eq. ([5]) that 
the nonlinear terms of both fields have opposite sign. 

It should be remarked that the results are obtained up to the second order 
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accuracy in Taylor expansion. In order to guarantee that the numerical solutions 
contain no large amount of truncation errors, the step sizes 5 and e are kept small 
enough. Nevertheless, the present method should still be good approximation to 
describe visually the mechanism of secondary protein folding. 

5. Conclusion 

An extension of phenomenological model describing the conformational dynam- 
ics of secondary proteins is proposed. The model is based on the matter interac- 
tions among relevant constituents, namely the nonlinear conformational field and 
the nonlinear sources. The fields are represented as the bosonic fields (/> and ip in 
the lagrangian. It has been shown that from the bosonic lagrangian with t/;^ self- 
interaction, the nonlinear and tension force terms appear naturally as expected, and 
coincide with some previous works ^ 10_ 

However the present model has different contour, and the folding process is 
getting slower since the EOMs governing the whole dynamics are the nonlinear Sine- 
Gordon and nonlinear Klein-Gordon equations. It is argued that the nonlinearity 
of the both fields are against each other. Note that the model is a generalization of 
earlier models which deployed both linear, or the linear and nonlinear equations. 
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